Condensation of classical nonlinear waves 
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We study the formation of a large-scale coherent structure (a condensate) in classical wave equa- 
tions by considering the defocusing nonlinear Schrodinger equation as a representative model. We 
formulate a thermodynamic description of the condensation process by using a wave turbulence the- 
ory with ultraviolet cut-off. In 3 dimensions the equilibrium state undergoes a phase transition for 
sufficiently low energy density, while no transition occurs in 2 dimensions, in analogy with standard 
Bose-Einstein condensation in quantum systems. Numerical simulations show that the thermody- 
namic limit is reached for systems with 16 3 computational modes and greater. On the basis of a 
modified wave turbulence theory, we show that the nonlinear interaction makes the transition to 
condensation subcritical. The theory is in quantitative agreement with the simulations. 

PACS numbers: 42.65.Sf,05.45.-a,47.27.-i 
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The problem of self-organization in conservative sys- 
tems has generated much interest in recent years. For 
infinite dimensional Hamiltonian systems like classical 
wave fields, the relationship between formal reversibil- 
ity and actual dynamics can be rather complex. Nev- 
ertheless, an important insight was obtained from soli- 
ton dynamics, in particular through numerical simula- 
tions performed for the focusing nonintegrable nonlin- 
ear Schrodinger (NLS) equation y], Q . These studies re- 
vealed that this system would generally evolve to a state 
containing a large-scale coherent localized structure, i.e., 
solitary-wave, immersed in a sea of small-scale turbulent 
fluctuations. The solitary-wave plays the role of a "sta- 
tistical attractor" for the system, while the small-scale 
fluctuations contain, in principle, all information neces- 
sary for time reversal. Importantly, the solitary-wave 
solution corresponds to the solution that minimizes the 
energy (Hamiltonian), so that the system tends to relax 
towards the state of lowest energy, while the small-scale 
fluctuations compensate for the difference between the 
total energy and that of the coherent structure @|. 

It is only recently that a statistical description of this 
self-organization process has been elaborated 0, Q . Re- 
markably it turns out that when the Hamiltonian sys- 
tem is constrained by an additional invariant (e.g., the 
mass), the increase of entropy of small-scale turbulent 
fluctuations requires the formation of coherent structures 
to "store" the invariant jj], so that it is thermodynami- 
cally advantageous for the system to approach the ground 
state which minimizes the energy Q. 

In the present work we consider the defocusing regime 
of the NLS dynamics, which is also relevant to the de- 
scription of thermal Bose gases [U El 0- This regime 
would be characterized by an irreversible evolution of 
the system towards a homogenous plane- wave 0, II , 
as described by the weak turbulence theory ^(J- This 
evolution is consistent with the general rule discussed 
above 0, IH 3 , because the plane- wave solution realizes 
the minimum of the energy (Hamiltonian) in the defocus- 



ing case. Accordingly, the NLS equation would exhibit 
a kind of condensation process, a feature that has been 
accurately confirmed only recently in the context of ther- 
mal Bose fields, where intensive numerical simulations of 
the NLS equation have been performed in 3D 0, Q • 

In this Letter we formulate a thermodynamic descrip- 
tion of this condensation process. We show that the 3D- 
NLS equation exhibits a subcritical condensation process 
whose thermodynamic properties are analogous to those 
of Bose-Einstein condensation in quantum systems, de- 
spite the fact that the considered wave system is com- 
pletely classical. Our analysis is based on the kinetic 
theory of the NLS equation, in which we introduce a 
frequency cut-off to regularize the ultraviolet catastro- 
phe inherent to ensembles of classical waves (Rayleigh- 
Jeans paradox). This allows us to clarify the essential 
role played by the number of modes in the dynamical 
formation of the condensate. Furthermore, we find that 
in 2D the NLS equation does not undergo condensation 
in the thermodynamic limit, in complete analogy with a 
uniform, ideal Bose gas. The significance of this result 
is that the system irreversibly evolves towards a state 
of equilibrium (maximum entropy) without generating a 
coherent structure that minimizes the energy (Hamilto- 
nian), in contrast with the general rule commonly ac- 
cepted HUSH. 

Given the universal character of the NLS equation in 
nonlinear wave systems, the reported condensation pro- 
cess could stimulate interesting new experiments in var- 
ious branches of physics. Nonlinear optics is a natural 
context where classical wave condensation may be ob- 
served and studied experimentally, in relation to the re- 
cent proposal of condensation for photons More- 
over, the formal reversibility of the condensation process 
could be demonstrated by means of an optical phase- 
conjugation experiment 12]. In addition to optics, wave 
condensation could be relevant for hydrodynamic surface 
waves on the basis of the recent progress on mea- 
surements of Zakharov's spectra. Moreover, our work is 
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also relevant to the important issue of strong wave turbu- 
lence 0| , where the emergence and persistence of large 
scale coherent structures in the midst of small scale fluc- 
tuations is a common feature of many turbulent fluids, 
plasmas and optical wave systems. 

We consider the normalized defocusing NLS equation 
in D spatial dimensions for the complex function t/>|15|: 



that this cut-off arises naturally in numerical simulations 
through the spatial discretization of the NLS equation, 
and manifests itself in real physical systems through vis- 
cosity or diffusion effects at the microscopic scale [lCj . 

To begin let us analyze the equilibrium distribution (3J 
in 3D. From energy and mass conservation one gets 



id t ip = -Aip + \ip\ 2 ip 
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where A stands for the Laplace operator in dimen- 
sion D. This equation describes the evolution of de- 
focusing interacting waves through the cubic nonlin- 
ear term. The dynamics conserves the mass (particle 
number) N = J \ip\ 2 d D x, and the total energy H = 

/(IWf + ihM 4 ) d D x. 

We address the dynamical formation of the conden- 
sate starting from a non-equilibrium stochastic initial 
condition for the wave ip, which we take to be of zero 
mean and statistically homogenous. In spite of the for- 
mal reversibility of the NLS Eq.Q), the nonlinear wave 
%p may exhibit an irreversible evolution towards ther- 
mal equilibrium. The salient properties of this evolution 
may be described by weak-turbulence theo ry. This the- 
ory derives, within certain approximations |l0| . an irre- 
versible kinetic equation for the averaged wave-spectrum, 
(flkiflks) = 0*1 — k- 2 )' ° k being the Fourier trans- 
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form of ip defined by a\^{t) — j ^(x, t) e 



— zk-: 



d t n kl (t) =Coll[n kl ] = Jd D k 2 d D k 3 d D k i Wk 1 ,k 2 ;k 3 ,k i 
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with W klM , HM = I ^ T 5W(k 1 +k 2 -k3-k 4 ) 5W(kl 

k 2 ~ k 2 — k 2 ) [ljj. As f° r the usual Boltzmann's equa- 
tion, Eq.J2J) conserves the mass N = V f rik(t)d D k, the 
kinetic energy E — V J k 2 n-^(t)d D k (being V the sys- 
tem volume for D = 3 and the surface for D = 2), 
and exhibits a iJ-theorem of entropy growth dS/dt > 0, 
where S(t) = J\n(n^)d D k is the non-equilibrium en- 
tropy. Accordingly, the kinetic equation © describes an 
irreversible evolution of the wave-spectrum towards the 
Rayleigh- Jeans equilibrium distribution |10| : 



These equations should be interpreted as follows. The 
initial non-equilibrium state of the field if) is characterized 
by a mass N and an energy E. The wave spectrum then 
relaxes to the equilibrium distribution (3J, whose tem- 
perature T and chemical potential /i are determined by 
Eas. l|4l5|) : a given pair (N,E) then determines a unique 
pair (T, /z). An inspection of Eq.Q) reveals that /z tends 
to 0, i.e. A c diverges to infinity, for a non-vanishing tem- 
perature T (keeping N/V constant), or a finite density 
N/V (keeping T constant). By analogy with the Bose- 
Einstein transition in quantum systems, this reveals the 
existence of a condensation process. The same conclusion 
follows from the analysis of the energy per particle E/N. 
There exists a non-vanishing critical energy per particle 
Etr/N = k 2 /3 such that n = 0. Dividing © by gj), one 
gets (E - E tI )/(Nk 2 ) = n^JI/(6k c ) + O (-fi/k 2 ). De- 
creasing the energy per particle effectively cools the sys- 
tem, and one reaches a finite threshold E tT , below which 

ni°¥tW) tion 0CCUrS ' 

"The same analysis in 2D readily gives N/V — 

+ 7rTln(l - k 2 /n) and E/N = fj, + k 2 /\n(l - k 2 /fi). In 
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where T and /x(< 0) are, by analogy with thermodynam- 
ics, the temperature and the chemical potential respec- 
tively. The spectrum @ is lorentzian and the correlation 
length of the wave ip is determined by /i, A c ~ l/s/—fi. 

The distribution © realizes the maximum of the en- 
tropy <S[nk] and vanishes exactly the collision term, 
Colling] — 0. However, it is important to note that 
Eq.Ql is only a formal solution, because it does not lead 
to converging expressions for the energy E and the mass 
N in the short- wavelength limit k — > oo. To regularize 
this unphysical divergence, we introduce an ultraviolet 
cut-off k c , i.e., we assume Uk(i) = for k > k c . Note 



this case /j, reaches : (i) for a vanishing temperature 
T (N/V constant), (ii) for a diverging density N/V (T 
constant), (Hi) for a vanishing energy per particle E/N. 
It results that condensation no longer take place in 2D. 
We have confirmed this result by direct numerical simu- 
lations of the wave equation (JJ. In contrast to the 3D 
case where the field ip generates a coherent plane-wave 
(condensate), in 2D ip remains a stochastic field of zero 
mean, whose spectrum evolves towards the equilibrium 
distribution (JjJ with fj, ^ Q, i.e., with a finite correlation 
length A c . 

The dynamical formation of the condensate has been 
studied using self-similar solutions of the kinetic equa- 
tion (5J, leading to a finite time singularity for some 
time £» 0, . Although this singularity is naturally 
regularized in the NLS dynamics, it describes explic- 
itly particle cumulation to zero wave- number k = 0. 
When the condensate begins to form (t > t*), an ex- 
change of mass between the condensate and the inco- 
herent wave-component is necessary to reach equilib- 
rium. It was argued in Refs. 0, |jj that this dynam- 
ics may be described by a three-wave kinetic interac- 
tion (the fourth wave being the condensate). This as- 
pect has been recently expressed more precisely by ex- 
tending the kinetic equation J2J to singular distributions 
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nt (t) = no(£)<5 (3) (k)+0 k (i) HE , so that a pair of cou- 
pled kinetic equations for the evolution of the condensate 
(no) and the incoherent wave-component (</>k) has been 
derived 01- These equations describe a flux of mass 
from the incoherent component towards the condensate, 
until thermal equilibrium is reached, i.e., until the colli- 
sion terms exactly vanish. This occurs for the equilibrium 
distribution (f)^ — T/k 2 , which actually corresponds to 
the distribution © with zero chemical potential. This 
allows us to legitimately assume p = below the transi- 
tion threshold E < E tr . 

The number of condensed particles no and the energy 
E(< En) may then be calculated by setting p, — in the 
equilibrium distribution One readily obtains (N — 
i2 )/V = AirTk c and E/V = A-KTk^/S, which gives 



n n /N 



n /N = 1 - E/Et, 



or n /N = 1 - T/T tI , where T tr = 3E tT /(4nVk^). As 
in standard Bose-Einstein condensation, no vanishes at 
the critical temperature T tr , and no becomes the total 
number of particles as T tends to 0. The linear behavior 
of no vs. E in Eq.JHJ) is consistent with the results of 
numerical simulations (see Fig. 1). However note that 
Eq.© is derived for a spherically symmetric continuous 
distribution of n k , while in the numerics the integration 
is discretized. Eq.© should thus be replaced by 
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N 
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N 



k y 
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(7) 
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where £\ denotes a discrete sum for —k c <k 
k c that excludes the origin k x = k y — k z = 0. 

This distribution is plotted in Fig^Jand compared with 
that given by the numerical simulations of Eq.Q. In 
our simulations we started with a non-equilibrium dis- 
tribution ip(x.,t = 0) = 27J fc akexp(zk • x), where the 
phases of the complex amplitudes a k are distributed ran- 
domly 0, Q ■ The numerical simulations confirm the ex- 
istence of the condensation process for sufficiently low 
energy densities 6]. We performed simulations with dif- 
ferent number of computational modes (8 3 , 16 3 , 32 3 , 64 3 
and 128 3 ). Our numerical results reveal that once the 
number of modes exceeds 16 3 , it only weakly affects the 
condensate fraction n^/N. This means that the system 
has reached some thermodynamic limit with only 16 3 
modes. Moreover, we verified that for a given total vol- 
ume V and a fixed density N/V, the energy threshold for 
condensation increases, as the frequency cut-off increases, 
in agreement with the theory (E tI oc k 2 ). 

The linear dependence Q between no and E gives a 
poor approximation of the numerical results, mainly be- 
cause the condensate fraction has been calculated by tak- 
ing into account only the linear kinetic contribution E to 
the total energy of the wave H . To include the nonlin- 
ear (interaction) contribution, we adapt the Bogoliubov's 
theory of a weakly interacting Bose gas [l^ to the clas- 
sical wave problem considered here. We start from the 
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FIG. 1: Condensate fraction uq/N vs. total energy density 
(H) /V, where (H) = E + Eo, Eq being the condensate en- 
ergy [see Eq.©]. Points (o) refer to numerical simulations of 
(6) the NLS Eq.JTJ with 64 3 modes for a fixed density N/V = 1 
(k c — 7r). The straight line (i) [(ii)] corresponds to the con- 
tinuous Eq. JHJ [discretized Eq.£J] approximation. Curve (Hi) 
refers to condensation in the presence of nonlinear interac- 
tions [from Eq.JSJ], which makes the transition to condensa- 
tion subcritical, as illustrated in the inset (with 1024 3 modes). 



total energy H of the nonlinear wave 

k ki,k2,k3,k4 

where <5k is the Kronecker delta symbol. The Hamilto- 
nian may be decomposed into four terms, H = H0 + H2 + 
H3 + H4, depending on how the zero mode, clq — <2k=o, 
and non-zero modes, ak^cb enter the expansion: 



H 
H 2 



1 

2V 
1 

E 



(|a | 4 + 2|a | 2 (iV-|ao| 2 )) 



(fc 2 4 


l«o| 2 ^ 
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Ha = 2V E ( 2a oa kl a k2 a k3 +c.c.)(5 kl +k 2 -k 3 



kl,k2,k3 



l ' 

Ha = — ^ a k 1 kj k3 k4^ki4k 3 -k 3 -k4 5 

fci ,fe 2 ,k 3 ,k4 

where ^~j' fc excludes the k = mode. The kinetic equa- 
tion requires the Hamiltonian to be diagonal in quadratic 
terms. To this end, we apply the Bogoliubov's transfor- 
mation for the canonical variables 6k = "kOk — fkilk) 
with the condition |iik| 2 — |«k| 2 = 1 that preserves the 
Poisson's bracket relation {ak,Ok} = * m the 6k's ba- 
sis. Imposing that the quadratic t erm H 2 is diago- 
nal in this basis, we find Uk = 

-k 2 



L 



1 



Li. with Lk 



h\ and Vk = 



2p + u B (k))/p and 
H2 = J2k' u B(k)blbu, where p = n /V = \a \ 2 /V, and 
w s(fc) — + 2p k 2 is the Bogoliubov's dispersion re- 
lation that takes into account the nonlinear interaction. 
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Let us emphasize that, in the 6k 's basis, the kinetic 
equations describing the coupled evolution of the conden- 
sate (no) and the incoherent wave-component (<£>k), are 
similar to those derived in Ref.|l7j. but the linear disper- 
sion relation uj(k) = k 2 is replaced by the Bogoliubov's 
expression ws(fc). The equilibrium distribution turns out 
to be ip e k q = T/uj B {k), with (WK>) = ^ 9<5 ( k - k '): so that 
the uncondensed mass (in the &k's basis) reads 

The total averaged energy (H) has contributions from 
Ho, H 2 and H A 0: (H) = E + Efc^K" = 
E + TJ2' k 1, where E = [iV 2 + (N — n ) 2 ] is the 
energy of the condensate. The temperature T may be 
substituted from this expression by using Eq.©, which 
gives a closed relation between (H) and no, 

{H) = Eo+ , (9) 

ELP 2 + po)A4(fc)] 

This expression is in quantitative agreement with the nu- 
merical simulations of the NLS Eq.JIJ, without any ad- 
justable parameter (see Fig. 1). Expression @ remark- 
ably reveals that the nonlinear interaction changes the 
nature of the transition to condensation, which becomes 
of the first order. This subcritical behavior is more pro- 
nounced when the number of modes increases. 

In summary, by considering the defocusing NLS equa- 
tion as a model, we showed that a classical nonlinear 



wave exhibits a subcritical condensation process in 31?, 
while no transition occurs in 2D. In spite of the for- 
mal reversibility of the NLS equation, the condensation 
process manifests itself by means of an irreversible evo- 
lution towards a homogenous plane- wave (condensate), 
with small-scale fluctuations superimposed (uncondensed 
particles) , which store the information necessary for time 
reversal. We formulate a thermodynamic description of 
the condensation process, whose properties are analogous 
to those of standard Bose-Einstein condensation in quan- 
tum systems. However, caution should be exercised when 
drawing conclusions about condensation in real bosonic 
systems, because the NLS equation only describes highly 
occupied modes satisfying nk ^> 1, so that it cannot de- 
scribe any transfer of mass between scarcely occupied 
high-energy modes and the condensate 0, B H U21 • Such 
a transfer is crucial at least at the beginning of the growth 
of the condensate that starts with zero mass. 



To conclude we would like to mention that the equilib- 
rium distribution |J3J also describes the velocity spectrum 
distribution in 2D fluid turbulence [2(j • An analogy be- 
tween 2D fluid turbulence and our results may be help- 
ful in the context of decaying turbulence at very high 
Reynolds number [2l|. In this case, following Onsager's 
argument |22| . the large scale vortices formed by the in- 
verse cascade would play the role of the condensate. 
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